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We discuss the statistical properties of index returns in a financial market just after a major 
market crash. The observed non-stationary behavior of index returns is characterized in terms 
of the exceedances over a given threshold. This characterization is analogous to the Omori law 
originally observed in geophysics. By performing numerical simulations and theoretical modelling, 
we show that the nonlinear behavior observed in real market crashes cannot be described by a 
GARCH(1,1) model. We also show that the time evolution of the Value at Risk observed just after 
a major crash is described by a power-law function lacking a typical scale. 



I. INTRODUCTION 

The time evolution of basic indicators of a financial 
market as, for example, the volatility of a market index is 
often showing a non-stationary pattern [Q. In fact, tests 
for whether volatility fluctuations of the size observed in 
empirical time series could be due to estimation error 
strongly reject the hypothesis of constant volatility |^|. 

The volatility time series, both historical and implied, 
experiences a dramatic increase at and immediately af- 
ter each financial crash ||. After the crash, the market 
volatility shows a stochastic evolution characterized by a 
slow decay. 

A non-stationary time evolution of volatility directly 
implies a non-stationary time evolution of the returns of 
a market index. The purpose of this paper is to show that 
the dynamics of market index returns just after a crash 
presents a statistical regularity with respect to the num- 
ber of times the absolute value of index return exceeds a 
given threshold value. 

The discovered statistical regularity is analogous to a 
statistical law discovered in geophysics more than a cen- 
tury ago and known today as the Omori law [Q. The 
Omori law states that the number of aftershock earth- 
quakes measured at time t after the main earthquake 
decays with < as a power-law function. The observation 
of this specific functional form in the stochastic relax- 
ation of the market statistical indicators to its typical 
values implies that the relaxation dynamics of a market 
index just after a financial crash is not characterized by 
a typical scale. The lack of a typical relaxation scale has 
direct implication on the estimation of risk measures as, 
for example, the Value at Risk of a stock portfolio ||. 

The structure of the paper is as follows. In Section 2 we 
report on the empirical properties of aftercrash market 



index time series. The discussion concerns the 1-minutc 
Standard and Poor's 500 index returns recorded after the 
Black Monday financial crash. Section 3 shows that the 
GARCH(1,1) model is unable to describe the empirical 
properties of the previous section. Section 4 discusses the 
estimation of the Value at Risk during the time period 
immediately after a financial crash. The paper ends with 
some concluding remarks. 



II. EMPIRICAL PROPERTIES OF 
AFTERCRASH SEQUENCES 

In this section, we characterize the time series of ultra 
high frequency returns after a major market crash. The 
variable investigated is the one- minute logarithm changes 
of a financial index r(t). A direct characterization of the 
time evolution of this random process is extremely dif- 
ficult in the time period after a market crash. This is 
due to the fact that the aftercrash period is highly non- 
stationary because the financial market needs some time 
to be back to a " normal" period. In order to characterize 
the aftercrash return time series we make use of a sta- 
tistically robust method. Specifically, we quantitatively 
characterize the time series of index returns by investi- 
gating the number of times \r(t)\ is exceeding a given 
threshold value |(| in the non-stationary time period. 

A similar approach is used in the investigation of the 
number n(t) per unit time of aftershock earthquakes 
above a given threshold measured at time t after the main 
earthquake. This quantity is well described in geophysics 
by the Omori law pj. The Omori law n(t) cx t~ p says 
that the number of aftershock earthquakes per unit time 
measured at time t after the main earthquake decays as a 
power law. In order to avoid divergence at t = 0, Omori 
law is often rewritten as 



n(t) = K(t + T)- p , 



(1) 



t 
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where K and r are two positive constants. An equiva- 
lent formulation of the Omori law, which is more suitable 
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for comparison with real data, can be obtained by inte- 
grating Eq. (1) between and t. In this way the cumu- 
lative number of aftershocks after the main earthquake 
observed until time t is 

N(t) = K[(t + r) 1 -P-T 1 -r>}/(l-p), (2) 

when p^l and N(t) = K ln(</r + 1) for p = 1. 

When the process is stationary the frequency of af- 
tershock n(t) is on average constant in time and there- 
fore the cumulative number N(t) increases linearly in 
time. We have tested that N(t) increases approximately 
linearly in a market period of approximately constant 
volatility such as, for example, the 1984 year. For inde- 
pendent identically distributed random time series it is 
possible to characterize n(t) in terms of an homogeneous 
Poisson process 0] . 

As an example of the behavior of N(t) empirically ob- 
served after a market crash, we investigate the index re- 
turns during the time period just after the Black Monday 
crash (19 October 1987) occurred at New York Stock Ex- 
change (NYSE). This crash was one of the worst crashes 
occurred in the entire history of NYSE. The Standard 
and Poor's 500 Index (S&P500) went down 20.4% that 
day. In our investigation, we select a 60 trading day af- 
tercrash time period ranging from 20 October 1987 to 14 
January 1988. For the selected time period, we inves- 
tigate the one-minute return time series of the S&P500 
Index. The unconditional one-minute volatility is equal 
to a = 4.91 x 10~ 4 . In Fig. 1 we show the cumula- 
tive number of events N(t) detected by considering all 
the occurrences observed when the the absolute value 
of index return exceeds a threshold value £ = 4a. We 
observe a nonlinear behavior in the entire period. Fig. 
1 also shows our best nonlinear fit performed with the 
functional form of Eq. (2). The agreement between em- 
pirical data and the functional form of the Omori law is 
pretty good. With a single time series it is not possible 
to have a good description of n(t) due to the discreteness 
of the data. In order to solve this problem we compute 
n a v(t) which is the running average of n(t) in a sliding 
window of 200 trading minutes. In the inset of Fig. 1 we 
show this quantity as a function of time in a bilogarith- 
mic plot. A power-law behavior is observed for almost 
two decades and the exponent of the best power-law fit 
is in agreement with the exponent p obtained by fitting 
N(t). 

We observe a similar behavior when we set the thresh- 
old value £ equal to 5a, 6a and 7a. The best fit of the 
exponent p slightly increases when £ increases and even- 
tually converges to a constant value. 

This paradigmatic behavior is not specific of the Black 
Monday crash of the S&P 500 index. In fact, we observe 
similar results also for a stock price index weighted by 
market capitalization for the time periods occurring af- 
ter the 27 October 1997 and the 31 August 1998 stock 
market crashes [|| . This index has been computed select- 
ing the 30 most capitalized stocks traded in the NYSE 
and by using the high-frequency data of the Trade and 
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t (trading day) 

FIG. 1: Cumulative number N(t) of the number of times 
\r(t)\ is exceeding the threshold £ = 4a during the 60 trad- 
ing days immediately after the Black Monday financial crash. 
The parameter a is the standard deviation of the process r(t) 
computed over the entire investigated period. The dashed line 
is the best fit of Eq. (2). The value of the exponent p is 0.85. 
In the inset, we show n av (t) which is a moving average of n(t) 
in a sliding window of 200 trading minutes. The dashed line 
is a power-law curve with an exponent equal to —0.85 and is 
shown as a guide for the eye 

Quote (TAQ) database issued by the NYSE. 

Our empirical results show that index return cannot be 
modeled in terms of independent identically distributed 
random process after a big market crash. In ref. j(| we 
introduce a simple stochastic model which is able to ex- 
plain the Omori law in financial time series. Specifically 
we assume that the stochastic variable r(t) is the product 
of a deterministic time dependent scale j(t) times a sta- 
tionary stochastic process r s (t) during the time period 
after a big crash. Under these assumptions, the num- 
ber of events of \r(t)\ larger than I observed at time t is 
proportional to 

n(t) oc 1 - F s , (3) 

where F s (r s ) is the cumulative distribution function of 
the random variable r s (t) . 

By assuming that the stationary return probability 
density function behaves asymptotically as a power-law 

fs{r s ) ~ -Jpr, (4) 
r s 

and that the scale of the stochastic process decays as a 
power-law 7(f) ~ t^ 13 , the number of events above thresh- 
old is power-law decaying as n(t) ~ (j(t)/£) a ~ l/t p . 
The exponent p is given by 

p = ap. (5) 
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FIG. 2: The dotted line is the mean value of the cumula- 
tive number N(t) of the number of times \r(t)\ is exceeding 
the threshold £ — 4<r over 10 4 simulations obtained with a 
GARCH(1,1) model with parameters estimated in a period of 
60 trading days immediately after the Black Monday financial 
crash. The dashed line is the mean of N(t) over 10 4 simula- 
tions obtained with a GARCH(1,1) model with parameters 
estimated in 6 nonoverlapping windows of 10 days each. As a 
reference we show also N(t) for the real data already shown 
in Fig. 1 (solid line). 



The previous relation links the exponent p governing the 
number of events exceeding a given threshold to the a 
exponent of the stationary power-law return cumulative 
distribution and to the (3 exponent of the power-law de- 
caying scale. 

The hypotheses of our model are consistent with re- 
cent empirical results. In fact, a return probability den- 
sity function characterized by a power-law asymptotic 
behavior has been observed in the price dynamics of 
several stocks || |9| and a power-law or power-law log- 
periodic decay of implied volatility has been observed in 
the S&P500 after the 1987 financial crash @. In |§ we 
have tested the hypotheses of our model by measuring 
independently the exponents a and [3 from real data and 
we have shown that the relation between exponents de- 
scribed by Eq. (5) is satisfied in all the three investigated 
crashes. 



III. COMPARISON WITH GARCH(1,1) MODEL 

Here we compare the empirical behavior of a market 
index after a major crash with the predictions of a sim- 
ple autoregressive model, specifically the GARCH(1,1) 
model (ll[| described by the equation 

of = a + a 1 r^_ 1 + /3ier£_ l5 (6) 

where r t is price return and at is the volatility. To this 
end we estimate the best value of the parameters of the 
GARCH(1,1) model. The estimation has been performed 



with the G@RCH 2.3 package Q. The best estima- 
tion of the parameters on the time series of one minute 
logarithmic price change during 60 trading days after 
Black Monday crash are «o = 2.87 10~ 8 , ax = 0.38 and 
Pi = 0.54. With these parameters we generate 10 4 sur- 
rogate time series according to Eq. (6) and we compute 
the average behavior of N(t). To mimic the dynamics af- 
ter the crash, we set as an initial condition a large value 
of return in each realization. Fig. 2 shows N(t) ver- 
sus time for real data and GARCH(1,1) model (dotted 
line) with a threshold I — 4cr (the same as in Fig. 1). 
The GARCH(1,1) time series converges to its stationary 
phase very quickly and it is unable to show a significant 
nonlinear behavior. 

In order to take into account the non-stationary be- 
havior of the return time series after a crash we have 
performed a different analysis. Specifically we divide the 
60 trading day period in 6 nonoverlapping time intervals 
of 10 trading days and we estimate the G ARCH (1,1) pa- 
rameters specific for each interval. We then generate 10 4 
GARCH(1,1) surrogate time series using the estimated 
parameters for each subinterval. The mean value of N(t) 
as a function of time obtained with this procedure is also 
shown in Fig. 2 as a dashed line. In each subinterval 
N(t) increases linearly for almost all the period. This is 
due to the fact that the non-stationary part is extremely 
small compared with the time scale of the figure. 

Below we theoretically describe the behavior of N(t) 
expected for a GARCH(1,1) time series. First of all we 
note that the expectation value of of conditioned to the 
value ctq is given by |ll| 

EtfWl) = + PxY + a \-^ 1+ /f (7) 

1 - («i + fix) 

This equation shows that, under the assumption of fi- 
nite unconditional variance ax + Px < 1, the mean value 
of the scale of the process decays exponentially to the 
unconditional value in a G ARCH (1,1) time series. The 
decaying time is given by t = — l/ln(ai + Pi)- With the 
parameters estimated we have r = 12 trading minutes in 
the case considered in Fig. 2. 

Since the variance of the GARCH(1,1) model depends 
only on cvo and on ax + fix, we perform a set of numeri- 
cal simulations of GARCH(1,1) time series with the same 
value of ao and by keeping constant ax+Px = 0.92, which 
are the values observed in the real data. Different time 
series of simulations are obtained by changing the value of 
a\ which controls the leptokurtosis of the unconditional 
density. Fig. 3 shows the average of n(t) over 10 5 realiza- 
tions of a GARCH(1,1) model with ax = 0.02, ax = 0.18 
and ax — 0.38. This last value corresponds to the one 
estimated from the real data. In the first two cases the 
unconditional kurtosis is finite, whereas for ax — 0.38 the 
fourth moment diverges. The value of the threshold is set 
to £ = 4cr in all cases. The unconditional value of n(t) is 
larger for larger values of ax- This is expected because 
the unconditional densities have the same variance and 
the tail are fatter for larger values of ax ■ For all values 
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FIG. 3: Mean of n{t), the number per unit of time that \r(t)\ is 
exceeding the threshold I = 4<r, for 10 5 numerical simulation 
of three different G ARCH (1,1) models with the same value of 
ao and of ai + P\ = 0.92. The values of ai are 0.02 (dotted 
line), 0.18 (solid line) and 0.38 (dashed line). This last value 
is the one estimated from the real time series of Black Monday 
aftercrash period. The gray line is the prediction based on the 
approximation of Gaussian density. 

of a>i, n(t) decays to the stationary value in less than a 
trading day. The decay is essentially exponential and a 
best fit with a function n(t) = a + b e~ ct gives a good 
agreement. 

Fig. 3 also shows that the larger is ot\ (and therefore 
the unconditional kurtosis), the faster is the decay of n(t) 
for very small values of t. This empirical behavior can be 
modelled and hereafter we show how it is possible to find 
an estimate of n(t) for small values of t. Let us consider 
that at time t = the crash occurs and the return is r . 
At time t = 1 the density is Gaussian with zero mean 
and variance given by o\ = ao + (a\ + I3i)tq. However 
at subsequent times the density is no more Gaussian. As 
a first approximation we can assume that the density is 
Gaussian with time-dependent variance described by (7). 
This is a very rough approximation, which is valid only 
when a.\ is very small. The gray line in Fig. 3 is n gau (t) = 
er f c(l / y/ 2at) which is the function obtained by using 
this approximation. For large values of t, n(t) tends to 
a constant value as e~*/ 2r . The agreement is very good 
in the case a\ — 0.02. However, for larger values of a.\ 
a better approximation is needed to explain the results 
of numerical simulations. A better approximation can be 
obtained by taking into account higher moments. Even 
if at time t = 1 the density is Gaussian, at t = 2 the 
density is no more Gaussian. Nevertheless we can find 
the exact expression for the expected value of the first 
four moments at t = 2 conditioned to the information at 
time t = 0. One finds 

E (r 2 )=Eo(r 3 2 )=Q, (8) 
Eo(i$) = a + (ax + frjo? = s 2 , (9) 
E (4)=6alai + E (rl). (10) 



TABLE I: Simulated and theoretical values of n(t = 2) in 
a GARCH(1,1,) model for different values of Qi and for the 
two regimes described in the text. Case (i) corresponds to 
r ( ) = 3.4 10 -3 and case (ii) corresponds to ro = 1 10 -3 . The 
other parameters are: a = 2.87 10" 8 , £ = 4a = 2.4 10~ 3 . 



Ql /3l 


case (i) 
n 3i m(2) nt eo (2) 


case (ii) 
n sim (2) n teo (2) 


0.02 0.90 
0.18 0.74 
0.38 0.54 
0.58 0.34 


4.41 10 _i 4.41 10" 1 
4.34 10" 1 4.31 10" 1 
4.11 lO" 1 3.95 10" 1 
3.73 10" 1 3.33 10" 1 


1.11 10~ 2 1.15 W~* 

1.26 10" 2 1.40 10" 2 
1.75 10~ 2 2.26 10~ 2 

2.27 10" 2 3.75 10" 2 



Under the assumption that the density is not very dif- 
ferent from Gaussian at time t = 2, we can approximate 
the density with a first order Edgeworth expansion (lj] 

where He^x) is the fourth Chebyshev-Hermite polyno- 
mial and K4 (t) is the fourth cumulant of n . By integrat- 
ing p(x) in Eq. (O) between I and infinity we can find an 
expression for n(t = 2). By performing the integration, 
we observe two different regimes: (i) when I 2 < 3(ao(l + 
on + 0i) + [on +/3i) 2 7q), n(t = 2) decreases when a\ in- 
creases; (ii) when I 2 > 3(a (l +a± + ("l + Pi) 2r o)j 
n(t = 2) increases when a\ increases. Table 1 shows two 
numerical examples corresponding to the two regimes. 
The results of Table 1 show that the agreement between 
simulated values and predicted values of n(t) is quite 
good especially for small values of a\. The fact that 
for large values of a\ the agreement is less precise is due 
to the significant contribution coming from the moments 
higher than the fourth that have been neglected in the 
Edgeworth expansion of Eq. (|i~l|). One could improve 
the forecast of n(t) by extending the Edgeworth expan- 
sion to higher order and by using the approximate density 
so obtained to calculate a better estimation of n(t). 

The technique of Edgeworth expansion can also be 
used to forecast the value of n(t) for value of t larger 
than 2. Our numerical simulations show that for moder- 
ate values of a\ and/or short forecast horizons the predic- 
tions work quite well. For larger value of a\ (especially 
when the fourth moment of the unconditional density is 
infinite) the prediction is valid only for small values of t. 

These theoretical and numerical observations confirm 
that the GARCH(1,1) is unable to model the power law 
decay of n{t) observed in real time series after a major 
market crash. 



IV. VALUE AT RISK AFTER A CRASH 

The empirical evidence and the theoretical model 
shown in the previous sections have direct relevance for 
risk management. One of the most widely used measure 
of risk is the so called Value at Risk (VaR) f|. The VaR 
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is the most probable worst expected loss at a given level 
of confidence over a given time horizon. Here we indi- 
cate the probability density for index return r over the 
time horizon r as ,f T {f)- The VaR Av a R associated to a 
certain probability of loss PvaR is defined implicitly by 




In stationary market conditions, the VaR is determined 
by the the density f T (r) . In Section 2 we have discussed a 
simple market model able to describe quantitatively the 
non-stationary evolution of index return time series after 
a market crash. In this model the shape of the return 
density remains constant but its scale changes in time. 
By using the relation r(t) = j(t)r s (t), it is direct to show 
that the (instantaneous) VaR after a market crash is not 
constant but changes in time as 

Ava R (t) = tWA&h, (13) 

where Ay R is the constant VaR obtained by using the 
stationary probability density function f s (r s ). We have 
seen Q that the occurrence of the Omori law after a 
market crash is consistent with a time evolution of the 
scale 7(f) = c\t~P + &z. Eq. (13) indicates that after a 
big market crash the VaR decreases in time as a power- 
law with exponent (3 (close to 0.2-^0.3 in the investigated 
crashes) and eventually relaxes to a constant value. 



V. CONCLUSIONS 



Our empirical observations show that the statistical 
properties of index return time series after a major fi- 
nancial crash are essentially different from the ones ob- 
served far from the crash. Other examples of statistical 
properties of market which are specific of the aftercrash 
period have been observed in the investigation of cross- 
sectional quantities computed for a set of stocks before, 
at and after financial crashes |l5|, |l6). The time period 
just after a crash is characterized by a relaxation to the 
typical market phase that can be modelled in terms of 
a power-law decay of the typical scale of index returns. 
We have shown that this observation rules out the possi- 
bility that the statistical properties of the time evolution 
of index returns can be efficiently modelled in terms of 
simple autoregressive models, such as the GARCH(1,1) 
model, in the non-stationary period after a crash. Our 
modelling of the aftercrash period also suggests that the 
Value at Risk of a financial portfolio measured just after 
a financial crash evolves in a non trivial way character- 
ized by a power-law evolution lacking a typical scale. In 
summary we conclude that the time period after a major 
market crash is characterized by statistical regularities 
which are specific to such a time period and not well de- 
scribed by models which presents a typical time scale in 
their stochastic evolution. The authors thank INFM ASI 
and MURST for financial support. 
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